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Introduction 


In  the  course  of  analyzing  data  under  contract  AF19628- 
81-C-0141,  an  attempt  was  made  to  fit  certain  exponential 
functions  to  number  density  data  obtained  from  axial  scat¬ 
tering  spectrometer  probes. 

These  data  were  collected  in  supercooled  clouds  assoc¬ 
iated  with  aircraft  icing.  It  was  noted  that  the  distribu¬ 
tion  of  the  number  density  spectra  in  the  size  range  of  2  to 
30  microns  diameter  seemed  to  be  describable  by  the  Weibull 
function  or  the  Khrgian-Mazin1  form  of  the  gamma  distribution. 
The  Weibull  function  has  been  used  in  time-to-failure  proba¬ 
bility  analysis.  These  functions  follow: 


Further  discussion  of  this  Weibull  function  can  be  found  in 
Bethea  et  al.2 


Gamma  Distribution  of  degree  2  (Khrgian-Mazin) 1 


x>0  a>0 

x<0 


(2) 


where  N  is  the  distribution  function 


The  method  of  least  squares  analysis  can  be  used  to  de¬ 
termine  values  of  a  and  b  in  these  functions.  It  is  noted 
that  this  effort  will  require  extensive  use  of  numerical 
analysis  to  solve  the  resulting  simultaneous  non-linear 
equations.  In  order  to  save  computation  costs  it  was  decided 
to  utilize  the  maximum  likelihood  method  (MLM)  which  has  the 
advantage  of  reducting  most  solutions  to  a  closed  form.  The 
maximum  likelihood  method  (MLM)  is  described  by  Breiman3. 


Weibull  Distribution  Function 


In  order  to  fit  a  function,  f (x) ,  to  a  data  distribu¬ 
tion  (y) ,  one  mathematical  method  which  can  be  used  is  the 
MLM.  This  method  requires  that,  for  the  entire  data  distri¬ 
bution,  the  sum  of  the  products  of  each  observed  value  and 
the  corresponding  value  of  the  function 


In  order  to  fit  this  function  to  (x^,y^)f  where, 

*  channel  i  diameter 
y^  *  channel  i  normalized  number  density 


"a"  and  "b"  must  be  solved.  Since  the  ln[f (x) ]  and  f (x) 
are  maximized  for  the  same  values  of  "a"  and  "b",  and  since 
f(x)  is  a  product  of  functions  of  "a"  and  "b" ,  it  is  easier 
to  solve  for  those  values  of  "a"  and  "b”  which  maximize  ln[f(x) ]. 


Therefore,  since  the  MLM  implies  that  the  partial  deri¬ 
vatives  equal  zero,  "a"  and  "b"  can  be  solved  for  by  the  fol¬ 
lowing  equations: 


jr  l  y.ln(Nabx.  e  )  =  0 
i=l  1 


and 


3b*  i  =  1  yiln<Nab*. 


b-1e'aXi) 


=  0 


Solving  equation  (1) : 


n 


Ta  1f1  y1ln<Habjt1 


.  ,  -ax. 

b-1  i 


e  ~)  *  0 


^  j  i  —ax. 

.l1  (yi}  3a  (ln(Nabxi  e  1 )) 


n 


E  (y^)  -jr-  (lnN+lna+lnb+(b-l) lnx. 
i=l  1  3a  1 


n  l  k 

1  yi  *  *i>  -  0 
i=l  1  a  1 


11  1  n  K 

s  yt  ^  yixI 

i=l  1  a  i-1  1  1 


n 


n 


A  .E  yi  ■  E  yixj 

a  i*l  1  i=l  1  j 


0 


-ax 


Solving  for  "a"  yields: 


h-  o* 


(5) 


n 

Z 

i=l 


(e+lnx . -ax^lnx . ) 

D  111 


0 


1  n  n  n  . 

t-  Z  y.  +  I  y.lnx.-a  Z  y.x?lnx.  =  0  (6) 

D  i=l  1  i=l  1  1  i=l  1  1  1 


Now,  substituting  for  "a"  (from  equation  (5))  in  equation  (6) 
and  realizing  that  !?  y.  jt  o,  yields: 

I*  JL 

i=l 


9 


b  n 


£  y . lnx. 


1  l 


£  V4 


Z  y.x.lnx. 
i=l  1  1  1 

n  K 

£  y .  x. 
i=l  Yl  1 


Since  b  cannot  be  solved  for  in  closed  form,  a  numerical 
method  must  be  used.  One  such  method  is  the  Newton- Raphson 
method  as  described  by  Scarborough4.  In  applying  this  method, 
let  equation  (7)  equal  g(b) .  Since  the  expression 


£  y . lnx . 
i-1  1  1 


yi 


is  a  constant  (K)  for  a  given  distribution,  and  independent  of 
the  choice  of  a,  b: 

n  b 
£  y.x.lnx. 

g(b)  =  £  +  K  -  - 

b  n  , 

£  Y-x* 

iii  1  1 


Now,  to  solve  for  the  root,  b,  such  that  g(b)  =  0,  an  initial 
guess  for  g  (bQ)  is  made..  Next,  solve  for  bn+1  where: 


b  .  ,  =  b 
n+L  n 


^  fbn> 

q'(bn) 


and  where 


bn+^-bn|  <  k  .  The  object  now  is  to  define  e  small  enough 


so  that  b  =  bn+^.  Since  the  Newton-Raphson  method  is 
repetitive  by  nature,  g'(bn)  can  be  approximated  without 
significant  loss  of  accuracy: 


g'(bn, 


lim 

E->0  E 


Therefore,  substituting  this  approximation  for  g'(b  )  in 
equation  8: 


g(b)E 


n 


Jn+1  bn  “  g(bn+E) -g(bn) 


Letting  E  =  .01(bn)  yields: 


or 


g(bn)  (.01)  (bn) 

bn+l  =  bn  "  g(b  +.01b  ) -g(b  ) 


n  '  n' 


n 


bn+l  "  bn  ” 


.01(bn) 

g(1.01bn) -g(bn) 
g(bn) 


n+1 


*  b  - 
n 


•01(bn) 

g(1.01bn) 

9(bn> 


-  1 


With  "a"  and  "b"  defined,  equations  (5)  and  ( 9)  respec 
tively,  only  the  independent  factor  N  remains  to  be  solved. 


Now,  it  can  be  shown  that  determining  the  value  of  N  by 
way  of  the  MLM,  as  was  done  for  "a"  and  "b"  in  equation  3  and 
4,  would  result  in  the  expression: 


Since  this  expression  is  undefined  for  N,  another  method  must 
be  utilized  to  determine  N.  Therefore,  applying  the  method 
of  least  squares  to  equation  1  and  solving  for  N  using  partial 
derivatives  yields: 
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3N 
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(yi“ 


.  .  -ax^ 

Nabx^”1e  l) 2 
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b  b 

n  K  T  u,  i 

E  2(y.-Nabx*r  e  x) (-abxTe  x)  =  0 
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e 
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-2yiabx^  Xe 
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i=l 


K  1 
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Factoring  out 


2ab  yields: 


n 

E 

i=l 


yixi"le 


n  .  ,  -ax^ 

N (ab)  E  (x?  e  x)  2, 
i=l  1 


and  solving  for  N  yields: 
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III.  The  Khrgian-Mazin  form  of  the  Gamma  Distribution  Function 

The  gamma  distribution  derivation  is  accomplished  in  a 
manner  similar  to  the  Weibull  function  by  fitting  f(x)  to 
(xi,yi)  by  way  of  the  MLM.  Thus,  given  that: 


(  Na3  „  2  ”axi 

x>0 

)T  xi  6 

1° 

x<0 

f(x> 


solving  for  "a"  via  partial  derivatives  yields: 


(2) 


n  n  j  -ax. 

3a  .E.  <yi1«<Nf“Xi2e  *)  *  0 
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(ri  15  (lnN+lnSI  +  InV-aitj) 
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S  Yi  C(^)  (^-)  -  x.)  -  0 
i=l  1  a3  2  1 
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i=l  Yi  («  “  Xi 


)  -  0 


3  " 

—  h  y . 

a  .  ,  1 1 
1=1 


di  Vi 
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Solving  for  "a"  yields: 
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IV.  Analysis 


The  Weibull  and  gamma  distribution  functions  were  applied 
to  axial  scattering  spectrometer  probe  data  obtained  from  26 
discrete  data  gathering  flight  passes  made  by  the  Air  Force 
Geophysics  Laboratory's  instrumental  C130-E  aircraft  during 
the  1979-1980  winter.  The  instrument  probe  was  designed  to 
collect  data  in  the  2-30  micron  range  at  2  micron  intervals. 

The  resultant  15  data  values  (channels)  for  each  flight  pass 
comprised  the  raw  data  sets.  However,  due  to  the  excessively 
noisy  data  signal  present  in  channel  1,  these  data  were  eli¬ 
minated  from  the  analyses.  Thus  data  channels  2-15  were  used 
as  input  to  the  Weibull  and  gamma  functions. 

In  order  to  demonstrate  the  ability  of  each  function  to 
approximate  the  number  density  distribution  of  each  flight 
pass.  Figures  1-3  are  presented  to  compare  the  actual  number 
density  distribution  for  channels  2-15  to  the  modeled  distri¬ 
butions.  In  addition,  the  pertinent  meteorological  data  (li¬ 
quid  water  content  and  temperature)  for  each  flight  pass  are 
indicated  on  each  figure. 

The  distributions  from  the  three  flight  passes  illustrated 
in  Figures  1-3  were  chosen  to  typify  the  range  of  liquid  water 
content  (LWC)  and  temperature  values  encountered  during  the  26 
flights.  As  the  figures  show,  both  functions  approximate  the 
number  density  distributions  for  channels  2-6  with  similar  ac¬ 
curacy,  but  the  Weibull  generally  fits  the  data  better  than 
the  gamma  function  for  channels  7-9.  However,  as  the  number 
density  data  levels  off  for  channels  11-15,  only  the  gamma 
function  consistently  models  this  behavior  in  an  acceptable 
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FIGURE 


CHANNEL  NUMoEk 


FIGURE 


FIGURE  4 
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manner.  The  Weibull* s  repeated  tendency  to  approach  zero  at 
a  rapid  rate  for  channels  13-15,  11-15,  and  9-15  respectively 
necessarily  eliminated  it  as  the  modeling  function  of  choice. 
(Note:  number  density  values  which  are  less  than  the  lower 
limit  of  the  Y  axis  (1.0E-1)  were  not  plotted.) 

The  impact  of  this  behavior  on  the  overall  performance 

of  the  Weibull  function  is  demonstrated  in  Figure  4,  which 

depicts  each  functions '  log  RMS  value  for  each  of  the  26 

flight  passes.  As  can  be  seen  from  the  figure,  23  of  the  26 

-3  -1 

(88%)  gamma  RMS  values  are  less  than  or  equal  to  2.0  cm  \i  ; 
whereas  only  seven  of  the  26  (27%)  Weibull  RMS  values  are  in 
this  category. 

This  dissimilarity  in  the  performance  of  the  two  func¬ 
tions  prompted  several  analyses  of  the  data  to  determine  if 
some  relationship  or  correlation  factor  could  be  discerned 
between  a  given  data  distribution  and/or  the  meteorological 
conditions  (LWC  and  temperature) ,  during  which  the  data  dis¬ 
tribution  was  obtained.  Among  these  analyses  was  the  compa¬ 
rison  of  the  ratio  of  the  number  of  LWC  values  above  the 
median  LWC  value  to  the  number  of  LWC  values  below  the  median 
for  each  function.  The  same  comparison  of  ratios  was  made 
using  temperature  as  the  variate,  with  neither  comparison 
evidencing  a  correlation  between  the  magnitude  of  LWC  or 
temperature  values  and  the  performance  of  the  Weibull  func¬ 
tion.  This  lack  of  correlation  between  flight  pass  variables 

and  the  Weibull  performance  held  true  for  log  RMS  values  both 

-3  -1 

above  and  below  the  2.0cm  m  delimiter. 

In  addition,  since  the  Weibull  performed  poorest  for 
data  channels  11-15  (usually  values  <  100),  an  attempt  was 


made  to  determine  if  the  Weibull  was  sensitive  to  a  particu¬ 
lar  range  of  data  values  less  than  100.  Despite  the  fact 
that  the  Weibull  attained  log  RMS  values  very  close  to  those 
of  the  gamma  function  for  seven  flight  passes  where  channels 
11-15  were  in  the  7-70  range  (see  Figure  4,  passes  17-23), 
there  were  an  equal  number  of  instances  where  the  RMS  values 
deteriorated  for  flight  passes  where  channels  11-15  were  in 
this  same  range.  (Note:  the  two  flight  passes  (2  and  14)  for 
which  the  Weibull  attained  the  worst  RMS  values  were  considered 
an  inadequate  sample  on  which  to  perform  analyses.  However, 
it  can  be  mentioned  that  the  gamma  function  also  attained  its 
worst  RMS  values  for  the  same  two  passes,  albeit  an  order  of 
magnitude  difference  exists  between  these  values  for  the  re¬ 
spective  functions.) 


Conclusions 


These  observations  and  analyses  indicate  that  for  this 
application,  the  gamma  distribution  function  is  a  more  accu¬ 
rate  and  consistent  modeling  method  with  which  to  approximate 
the  data  distributions.  Subsequence  analysis  is  planned  to 
develop  a  weighting  scheme  to  refine  these  modeling  techni¬ 
ques  . 


References 


1.  Khrgian,  A.  Kh.,  Mazin,  I.P.,  Analysis  of  Methods  for 
Describing  Cloud  Droplet  Distribution  Spectra,  Trudy 
TsAO,  No.  17,  1956. 


2.  Bethea,  R.  M. ,  Duran,  B.  S.,  Boullion,  T.  L.,  Statisti¬ 
cal  Methods  for  Engineers  and  Scientists,  Marcel  Dekker 
Inc.,  1975,  p.  66. 


3.  Breiman,  L. ,  Statistics:  With  a  View  Toward  Applications, 
Houghton  Mifflin,  1973,  p.  194  (6.12). 


4.  Scarborough,  J.  B.,  Numerical  Mathematical  Analysis, 
Johns  Hopkins  Press,  1962,  p.  199. 


